Rotation Rate of Rods in Turbulent Fluid Flow 
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The rotational dynamics of anisotropic particles adverted in a turbulent fluid flow are important in many 
industrial and natural setting. Particle rotations are controlled by small scale properties of turbulence that are 
nearly universal, and so provide a rich system where experiments can be directly compared with theory and 
simulations. Here we report the first three-dimensional experimental measurements of the orientation dynamics 
of rod-like particles as they are adverted in a turbulent fluid flow. We also present numerical simulations that 
show good agreement with the experiments and allow extension to a wide range of particle shapes. Anisotropic 
tracer particles preferentially sample the flow since their orientations become correlated with the velocity gra- 
dient tensor. The rotation rate is heavily influenced by this preferential alignment, and the alignment depends 
strongly on particle shape. 



The dynamics of anisotropic particles in turbulent fluid 
flows are central to understanding many applications ranging 
from cellulose fibers used in paper making! to ice crystals in 
clouds 2 ' 3 - to locomotion of many micro-organisms 4 " 7 . When 
the particles are small and their concentration is low, their 
rotations are determined by the velocity gradient along their 
trajectories. Velocity gradients are dominated by the small- 
est scales in turbulence about which we have an extensive 
fundamental understanding^^. So the statistics of rotations 
of anisotropic particles forms an important problem founda- 
tional to many applications for which we need to develop a 
predictive understanding based on the fundamental properties 
of turbulence. 

Advances in imaging technology have made it possible to 
obtain time -resolved trajectories of particles in turbulence us- 
ing high speed stereoscopic imaging, and these measurements 
have produced many new insights about Lagrangian transla- 
tional dynamics of spherical particles 11 . There is also an ex- 
tensive literature on the motion of anisotropic particles in fluid 
flows. Anisotropic particles have fascinating dynamics even 
in simple flow*£iI£. Several simulations have addressed the 
turbulent casei^— , and experiments have measured dynam- 
ics in 2D flows 2 ! and orientation distributions in 3D turbulent 
flowsa 2 ^!. However, there are no available experimental mea- 
surements of time-resolved rotational dynamics of anisotropic 
particles in 3D turbulence. 

When ellipsoidal particles are small compared with the 
smallest length scales in the flow, their rotation rate is deter- 
mined by the velocity gradient tensor through Jeffery's equa- 
tion^: 

a 2 - 1 

pi = Qijpj + ^— (SijPj - PiPkSkiPi) (1) 

where pi is a component of the orientation director and a = 
l/d, is the aspect ratio of the ellipsoid given by the ratio of 
length (I) to diameter (d). f2,-j is the rate-of-rotation ten- 



sor, and Sij is the rate-of-strain tensor which are the anti- 
symmetric and symmetric parts of the velocity gradient tensor 
respectively. Particle dynamics become much more compli- 
cated when the particles are large-^ or density mismatched 2 , 
or high particle concentration produces inter-particle interac- 
tions and two-way coupling!. A phenomenology of the dy- 
namics of anisotropic tracers at low concentrations in turbu- 
lence is needed to provide a foundation for continued work on 
the more complex cases that often appear in applications. 

We have performed an experimental and computational 
study of rotation rates along the trajectories of anisotropic par- 
ticles in turbulent flows. In the experiments, small rods are 
tracked using stereoscopic images from four high speed cam- 
eras. Turbulence was generated between two grids oscillating 
in phase in an octagonal Plexiglas tank 2 ! that is 1 x 1x1.5 
m 3 . The experiments are performed at two Taylor Reynolds 
number, R x = {IbuL/v) 1 / 2 , R x =160 and R x =214. 
The Kolmogorov length scale (77 = (V 3 /(e}) 1/4 ) and time 
scale (t j( = (iV( e )) ) are respectively 77 = 375/im and 
T n = 70ms, and 77 = 210/itn and t v = 25ms. The rods are 
1 mm in length by 200/im in diameter and are cut from ny- 
lon thread at density of p = 1.15 g/cm 3 . The fluid is density 
matched with rods by adding 19% CaCl2 by mass to water. 
The position of the center of the rods was measured with an 
uncertainty of « 40/zm. We have measured the orientation 
of rods in two dimensional (2D) images from each camera 
and reconstructed the orientation of rods in three dimensional 
(3D) space using images from multiple cameras. Determin- 
ing the orientation of the rods in 3D requires measurements 
of the orientation in 2D images from at least three cameras 
that are not in the same plane. The rotation rate, p, is deter- 
mined from quadratic fits to the measured orientation versus 
time data. We have also studied the motion of anisotropic trac- 
ers in direct numerical simulations (DNS) of homogeneous 
and isotropic turbulence at R x =180. The translational mo- 
tion of tracer rods matches that of fluid particles, so we are 
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FIG. 1. Three-dimensional view of a rod trajectory with a large 
rotation rate from the experiment at R\ = 214. The color of the 
rod represents the rotation rate. This rod is tracked for 284 ms. The 
green lines show the projection of the center of the rod onto the y-z 
and x-y planes. The rod is a circular cylinder with length 1 mm and 
diameter 0.2 mm. 



able to use a database of previously simulated Lagrangian tra- 
jectories 2 * to integrate Jeffery's equation ([l} and obtain the 
time evolution of particle orientations. The spatial spectral 
resolution was 512 3 points. These simulations stored the full 
velocity gradient tensor along Lagrangian trajectories at time 
intervals of about l/10r r; . Equation Q] has been therefore in- 
tegrated a posteriori with the same integration time step and 
with an Adams-Bashforth second order in time scheme. 

Figure Q] shows an experimentally measured trajectory of a 
1mm rod at R\ = 214. This example illustrates several of 
the important properties of the rotation of rods. First, this rod 
has bursts of high rotation rate where the rotation rate squared 
is up to 30 times its mean reflecting the intermittency of rod 
rotations. Second, in the upper right, the rod is caught in a 
vortex, but its rotation rate is not large because the rod has 
become aligned with the vorticity reflecting the tendency of 
anisotropic particles to become aligned by the velocity gradi- 
ents in the flow. 

The probability distribution function (pdf) of the rotation 
rate squared, pip\, of rods is shown in Fig. [2^. The pdf has 
long tails which indicates the presence of rare events with ro- 
tation rates squared up to 60 times the average value ((piPi)). 
The agreement between DNS and experiment is very good for 
the core of the distribution up to 20 (piPi). At larger rotations 
rates the experimental pdfs are slightly below the DNS. This 
difference is not much larger than the systematic errors in the 
experimental data represented by the error bars, but it may re- 
flect the effect of finite length of the rods in the experiment. 
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FIG. 2. The pdf of rotation rate squared of rods, (a) Comparison of 
simulation (Blue line, R\ — 180) with experiment (Black o, R x = 
160; Red +, R\ — 214) for aspect ratio a = 5. Error bars represent 
both random uncertainty and the systematic uncertainty produced by 
measuring rotation rates from experimental orientation data over a 
time interval, (b) Comparison of the pdf of rotation rate squared for 
rods at a = 5 (Blue +) and spheres at a = 1 (black o) with pdfs 
of enstrophy (Green □) and energy dissipation rate (red >) for the 
simulation at R\ — 180. One symbol is displayed for every twenty 
bins. 



The rod lengths are 2.677 at i? A =160 and 4.7r] at R X =2U. Ref 
[16] indicates that rods less than about 7?/ should not have 
a measurable change in their rotation rate variance from the 
tracer limit, but it is possible that the rare events are more sen- 
sitive to rod length. In Fig.|2j), the probability distribution of 
rotation rate squared from the numerical simulation is com- 
pared with the distributions of enstrophy (uj 2 = 2f2yf2,j) and 
energy dissipation rate (e = 2vSijSij, where v is the kine- 
matic viscosity). For large values of rotation rate squared, the 
pdf for rods (a = 5) is larger than that of spheres (a = 1), and 
both are larger than either enstrophy or energy dissipation rate. 
Spheres simply rotate with the vorticity, pipi = fiy flij sin 2 9, 
and the larger intermittency for spheres compared with enstro- 
phy comes from the distribution of the angle, 8, between p and 
the vorticity vector. Rods with a = 5 have a contribution to 
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their rotation rate from Sij, which one might think would de- 
crease the probability of large rotations, making it more like 
the distribution of energy dissipation rate. However, varia- 
tions in the orientation of the rod with respect to the velocity 
gradient tensor contribute additional fluctuations giving rods 
the most high rotation rate events. 

Figure [3] shows the effect of the shape of particles on their 
rotation rate in turbulence. The rotation rate variance of disk 
shaped particles (a < 1) is much larger than that of the 
spheres (a — 1). This can be qualitatively understood as the 
additional contribution of strain (Sij in equation [TJ to the ro- 
tation rate. However, the rotation rate of rods (a > 1) is much 
smaller than spheres even though the rate-of-strain contributes 
to their rotation as well. At large a, our simulations agree with 
earlier work on thin rods at lower Reynolds number by Shin 
and Kochi£. The experimental measurements at a = 5 are 
consistent with the simulations considering the measurement 
uncertainties. 

Understanding the rotation rate data in Fig. |3]requires con- 
sidering the preferential alignment that occurs between parti- 
cles and the velocity gradient tensor. When particles are ori- 
ented randomly, their rotation rate variance can be calculated 
analytically from Eq. Q] by extending the calculation in Ref 
[16] to finite aspect ratio: 
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This result is shown as the green solid curve in Fig [3^. As par- 
ticles are advected by the flow, they become oriented so that 
their rotation rates are very different than the randomly ori- 
ented case, with the largest difference occurring for thin rods 
(a >> 1). Thin tracer rods are material line segments, so the 
large decrease in rotation rate for thin rods can be qualitatively 
understood as resulting from the known preferentially align- 
ment of material lines with the vorticity vecto r 18 ' 27,28 . Since 
only the vorticity perpendicular to the rod contributes to its 
rotation rate, aligned particles have greatly reduced rotation 
rates. 

Figure [3j? shows the dependence of rotation rate on aspect 
ratio for several simple cases in order to better understand the 
ways the data in Figure [3^ depends on both alignment with 
the velocity gradient tensor and the time evolution of the ve- 
locity gradient tensor. If the velocity gradients are sampled 
from DNS but are delta-correlated in time, the particles do not 
develop alignment with the velocity gradients, and therefore 
the rotation rate variance matches the prediction for randomly 
oriented rods in Eq. [2] When the velocity gradients are sam- 
pled from the DNS and held fixed in time, the alignment effect 
is too strong and the rotation rates are smaller than the rotation 
rates of rods advected in the turbulence at almost all aspect ra- 
tios. Finally, we integrated particles through time independent 
plane shear flow which generates Jeffery orbits. None of these 
cases has even qualitative agreement with the rotation rates of 
rods advected in the turbulence. Since the dependence of the 
rotation rate variance on aspect ratio is sensitive to both in- 
stantaneous statistics and temporal correlations of the velocity 
gradient tensor, it provides a metric for evaluating models of 



the velocity gradients in turbulence that is directly accessible 
experimentally. 
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FIG. 3. Rotation rate variance as a function of aspect ratio, (a) Blue 
+ is the DNS at R x = 180. Black □ is the experiment at R x = 214. 
Red o is the experiment at Rx — 160. The error bars on the exper- 
imental points represent systematic error due to extrapolation from 
the finite fit time required to measure the rotation rate~ The pur- 
ple o shows the result for infinite aspect ratio rods from simulations 
by Shin and Kocb~ at Rx = 53.3. The Green line is the analytic 
prediction for randomly oriented rods from Eq. |2] The dashed line 
indicates the rotation due to vorticity. (b) Simple cases of particles 
integrated through different velocity gradient fields: Blue + uses ve- 
locity gradients along Lagrangian trajectories at Rx = 180. Black □ 
uses velocity gradients sampled from the DNS and delta-correlated 
in time. This data matches the analytic prediction for randomly ori- 
ented rods shown as the green line. Red • uses velocity gradients 
sampled from the DNS and held fixed in time. The purple -4 uses the 
velocity gradient of a plane shear flow, (c) Comparison of simula- 
tions at different Reynolds numbers. 



Figure [3J: shows the dependence of the rotation rate on the 
Reynolds number. In simulations ranging from Rx = 31 to 
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FIG. 4. Fourth moment of the rotation rate of ellipsoids. Blue + is 
the DNS at R x = 180. Black ■ is the experiment at R\ = 214. 
Red o is the experiment at R\ — 160. Green dashed line is the 
prediction for spheres uncorrelated with the direction of fiy. The 
purple dashed-dot line is the fourth moment of Sij and the black 
dashed line represents the case of a Gaussian distribution. 



R\ = 400, the rotation rate variance changes by less than 9% 
with the largest differences occurring for disks (a << 1). 
For a more quantitative evaluation of the shape of the pdf 
of rotation rate squared, we report in Fig. a normalized 
fourth moment of the rotation rate, ((piPi) 2 ) / (piPi) 2 as a 
function of aspect ratio. The experimental measurements are 
in fairly good agreement with the simulations considering that 
the differences are on the order of the experimental measure- 



ment errors and particle size may affect the tails of the ex- 
perimental distribution as discussed above. The fourth mo- 
ment for a sphere can be related to the fourth moment of 
the vorticity tensor by assuming that the orientation director 
is uncorrelated with fiy which leads to ({piPi) 2 ) / (piPi) 2 = 
6/5((fiyfiy) 2 )/ (ClijQij) 2 , in good agreement with the sim- 
ulations. All these fourth moments are much larger than the 
value of 5/3 that is obtained if the components, p\, have a 
gaussian distribution. Rods (a >> 1) and disks (a << 1) 
have nearly identical normalized fourth moments, and the 
variation with aspect ratio is less than 20%, indicating that 
the dependence on particle shape of the normalized probabil- 
ity distribution is much weaker than that of the variance of the 
rotation rate in Fig [3] 

From both experiments and numerical simulations we have 
obtained a phenomenology of the rotational dynamics of 
anisotropic particles in turbulent fluid flow. Rod-like particles 
have a rotation rate that is strongly affected by the alignment 
of the rods with the vorticity vector. Disk-like particle also 
show effects of alignment although their rotation rate vari- 
ance is closer to the randomly oriented case. The transition 
between the two limiting aspect ratios is quite abrupt, with 
80% of the difference in rotation rate variance between disks 
and rods occurring between aspect ratios 0.5 and 2. Thus for 
many bacteria or ice crystals that are transported in turbulent 
flows, the picture of a thin Lagrangian rod or disk is much 
better than the approximation that they are spheres. 
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